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Abstract 

We introduce new techniques that can preserve unitarity of the system including ghost 
particles. Negative norms of the particles can be involved in zero- norm states by 
constraints of the physical space. These are useful to apply the higher-derivative 
propagator for quantum gravity to suppress divergences of vacuum energy and graviton 
mass correction. The quantum effects are mainly depending on the ghost mass scale. 
As the scale can be chosen in any order, the observed cosmological constant is realized. 
Further, applying ghost partners for the standard model particles, quantum gravity 
with matter fields becomes renormalizable with power counting arguments. 



1 Introduction 

There is a long way to unify general relativity and quantum mechanics as in the framework 
of quantum gravity [TJ EJ |3] . The two leading candidates of quantum gravity are string 
theory and loop quantum gravity. Unfortunately, as any theory of quantum gravity has 
some deep problem [3], the theory is not complete yet. In general, it is usually thought 
that something new must happen at the Planck scale to make a consistent theory of 
quantum gravity. Nevertheless, the scale is too high to reach with current experiments, 
hence, there is no experimental hint for the statement. Meanwhile, considering the cutoff 
scheme, observed cosmological constant suggests that the cutoff scale of vacuum energy 
should be around the neutrino mass or micrometer scale. To accommodate with the 
observation, new method that can have small energy scale seems necessary. 

Perturbative approach of quantum gravity is appropriate for low energy U |7]. 
Taking the flat-space background, one can quantize the weak gravitational field. In this 
method, there appear bad divergences in many Feynman diagrams because the coupling 
constant has negative dimension. For instance, quantum correction to the graviton mass 
squared has quartic divergence. If the cutoff scale is the Planck mass, graviton mass is 
correspondent to the Planck mass. This is similar to the Higgs mass problem though it 
is renormalizable and less problematic for physics. The problem of graviton mass is more 
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serious because observations of galaxy and clusters ensure a very long-range force but 
there is no natural way to protect the mass from quantum correction. 

Quartic divergence is the worst in the perturbative quantum gravity for one-loop. 
This is the same to the divergence of vacuum energy. Then, if some method can solve 
the cosmological constant problem, the problem of perturbative quantum gravity may 
be also solved. The simplest way is to take into account a ghost particle [5]. For the 
sake of opposite sign of commutation relation, divergences induced by normal particle 
(meaning it has positive norm) can be canceled. A difficult problem relating to the ghost 
is the violation of unitarity that indicates the violation of probability conservation due to 
negative norm. Although there are several discussions and possible solutions e.g. [9l [TO], 
the problem still remains in general understanding. This paper suggests new solution to 
avoid unitarity violation by selecting proper physical state relating between normal and 
ghost particles. In particular, the mixing state of them has zero norm that plays a key 
role to hold positive semi-definite norm. All the conditions of unitarity can be satisfied 
by the constraints to be physical states. 

We propose the existence of ghost partner for graviton. To be consistent with large- 
scale gravity, the ghost should have a mass. In order to realize the higher-derivative 
propagator, usual graviton is ought to be massive because there is a mismatch between 
the propagators of massless and massive gravitons. A mass term of graviton induces 
negative norm in general so that we choose the Pauli-Fierz mass term. Taking the term 
for both gravitons, all divergences can be improved and the mass squared correction 
of graviton will be better. However, to achieve consistent theory without fine-tuning, 
additional gravitons and their partners are necessary to increase the power of propagator. 
They are assumed to have the same quantum number of the original graviton but with 
different masses. If some conditions are satisfied, the theory can be super-renormalizable 
or completely finite. This can be seen as a modified version of the Lee- Wick model |12j . 

When we take enough number of gravitons and some conditions, quantum effects 
including vacuum energy are sufficiently suppressed. Amounts of the effects are mainly 
dependent on masses of gravitons instead of the cutoff scale. As the unitarity violation 
can be avoided, scales of the masses are not need to be very high such as the Planck 
mass. Hence we possess the correct value of the cosmological constant by taking the 
graviton masses around the neutrino mass scale. This implies that gravity is modified at 
micrometer scale. For the Yukawa-type modification, the inverse square law of gravity is 
confirmed on the scale longer than 56/mi at 95% precision [13] and more stringent bound 
about 20/im is recently given by [H] . It is not easy to be consistent with these experiments. 
On the other hand, the loop calculation yields small graviton mass which is not negligible 
even for the case getting the correct vacuum energy. Writing the cosmological constant A 
and the gravitational constant G, the lightest graviton mass becomes about y/GX which 
is of order 10~ 10 pc _1 . This is also close to the experimental bound of the lightest graviton 
mass [15]. Thus, the theory of quantum gravity without fine-tuning makes modifications 
for short and long distances simultaneously and both are near to current experimental 
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bounds. In the sense of possible experimental test, quantum gravity with multi-gravitons 
in low energy is interesting. 

The paper is organized as follows. In section 2, we explore the model of multi-gravitons 
with Pauli-Fierz mass terms, afterward, the gravitational wave for each graviton is canon- 
ically quantized. Section 3 introduces the interaction with fermions using connection and 
covariant derivative to calculate the gravitational potential. Section 4 explains how to 
preserve unitarity from ghost particles of negative norm by selecting physical state and 
wave packets. In section 5, vacuum energy is calculated with higher-derivative propagator 
mentioning the modification of gravity for short distance to keep it viable with experi- 
ments of gravity and the cosmological constant. In section 6, loop correction of lightest 
graviton mass is calculated with one contribution of a Feynman diagram. At last, we have 
a short conclusion at section 7. 

In our notation, we use if" 11 = diag(— 1, 1,1,1) and natural units c = h = 1. 



2 Quantization of gravitons and ghost gravitons 

Throughout this paper, we assume gravitational wave is quantized [16] which can be 
done in the same way of other fundamental fields, especially similar to photon. Though 
this is not a standard method to quantize gravity, the scheme is convenient to apply 
higher-derivative propagator and to see unitarity of the system. Some problems may arise 
since the quantization violates Lorentz invariance and leads unusual interaction terms for 
momentum operator. Further, we need to cancel out many terms inside the Lagrangian 
to quantize the gravitons. It looks like very schematic but can be a powerful method. As 
an effective field theory, the system will be valid. 



2 . 1 Lagrangian 

Let us start with the Einstein-Hilbelt action 

£ = + £mat, (1) 

167rG 

which leads the Einstein equation G^ v = RnGT^. In the Minkowski space background, 
the metric can be expanded as 

= rip, + «V> 9^ = rT - ^ + K 2 h%h» v . (2) 
Similarly, yj—g can be written as 

=l + - K h- -K 2 h^ u h% + \n 2 hh + • • • , (3) 
2 4 o 

where h = r\ ilv h iXV . To make the higher-derivative propagator, we explore multi-gravitons 
model. The gravitons are denoted by h$), where n runs from zero to M — 1. These fields 
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are dealt with a part of the real gravitational field i.e. h pu = h$ . All the gravitons 
are assumed to be massive with Pauli-Fierz mass terms 

Anas = - E(-l) n <(^> (nW - h^h^) (4) 



n=0 



Mass terms of ghost gravitons take opposite signs, and relating to those minus signs, ghost 
kinetic terms are modified by 



.2 



Af/2-1 



j2 (h( 2n+1 ^nh% n+ v - 2h% n+1 ^d"d a h^ n+1 ^ 



K 

~ \MGJ=g~ ^ v " ""^ ""^ " (5) 

v a n=0 v ' 

+ 2h^d^d u h^ 2n+1 ^ - h^+vnh^+V)), 

where we ignore total derivatives. Still, it is not enough to make canonical quantization 
for each field h$ . To remain only necessary terms, following terms are taken for the 
cancellation: 

1 K 2 * 

£ kin = - -= V (/i (ra) ""D^; ) - 2h {m) ^d'"d (T h {n)au 

J—g32n& F ia\ 

+ 2h {m) ^d ll d v h {n) - h {m) Uh {n) )). 
The energy momentum tensor is given by = -^{d p ^^^y- — ^g^— then 

(-l) n 8nGT^ J-K(d»d,h^ + &>d v h„M - - d,d u h^) 

1 1 (7) 

- -K(d p d a h^ - Uhfri)^ + ~Km 2 (h<$ - h^ Vllu ). 

Multiplying for both sides and applying d^T pv = as the first order approximation, 
we get d^h$ — d u h^ = 0. Plugging them into the equations above, we get 

(□ " m*)h$ = -(-l)^(T^ + _ V , U )T). (8) 

They can be related to the equations of motion of higher- derivative theory. Considering 
Af = 2, they yield 

(□ - ml)(D - m\) V = K - mg) — (T^ + ^(^4 - ^)T). (9) 
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2.2 Quantization 

Hereafter, we choose the transverse-traceless gauge (TT-gauge) for all gravitational fields. 
If gravitational wave of each field can be quantized we can write 




where E p = \/p 2 + ni\. Creation and annihilation operators obey 

a^\0) = 0, [a^),a^' X)t ] = {-l) n 5 nn ,5 xx ,8\p - p'). (11) 
Polarization tensors have following properties [21] : 

$ = 42, e (A)M M = 0, p^ = 0, e We^ = ^'. (12) 
In the TT-gauge, the Einstein tensor with second order of expansion can be calculated as 

Gf u = -K 2 (h pa d,d„h^ + h w uh" v + x -h pa nw° %w ). (13) 

Note, the last two terms are automatically canceled if graviton is massless. To perform 
the canonical commutation relations for massive gravitons, they are modified by taking 

c can = —-^=( ^'s^^w^ + ^ v^ D V + ^V D ^)- ( 14 ) 

As a convention, we set VlQnG. The energy momentum tensor now becomes T pu = 
(d^hpa d u h^ pa ) . For the quantization of gravitational wave, the Hamiltonian is given 
by H = J d 3 x T 00 then 

where we omit the zero point energy. This is a familiar form and commutation rela- 
tions are as usual except minus sings: [H,a p X ' n ^\ = {—l) n E p n ^a p X ' n ^. and [P, ap A ' n ^] = 
(-l) n pa { p ' n) \ where P = J d 3 xT 0i . Writing = (H,P), they hold the identities 

e -iP.x a (A,n)t ei p. x = a (\,n)i e - ( -i rip . x ^ e _ lP . x ^ ^ = ^.^-(-i)"^* ^ 

provide the consistent specetime dependence of h^)(x). In addition, energies of ghost 
particles are positive [T7] . 

Here, let us note the aspect of this quantization scheme. At first, the interaction part 
of the Hamiltonian violates Lorentz and CPT symmetries. In particular, matter-gravity 
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couplings are important to see the effect of the violations. Various scenarios are sought 
to observe the Lorentz violation [IB]. Other perturbative methods of quantum gravity 
such as the Lagrange formalism [5] and the ADM action [19] rarely violate the Lorentz 
symmetry so that it is a characteristic feature. Secondly, interaction terms appear for 
the momentum operator since G° l and T° l of matter sector in a curved spacetime are not 
zero in general. Adding to the time evolution, the configuration of the wave packet gains 
space evolution by the interaction via gravitons. It is not discussed in the literature and 
we do not know how it affects to particle physics. This effect may be observable as the 
size of wave packet is relating to the probability density. At any rate, the gravitational 
interaction is very weak so that it will not leave problematic result. Other important thing 
to discuss is how to derive mass, ghost kinetic, and cancellation terms. There needs more 
fundamental theory to produce them. Considering the theory including higher curvature 
terms, they will appear naturally and the quantization can be performed without adding 
terms by hand. However, higher curvature terms make stronger divergences and we cannot 
avoid fine-tuning to cancel out them. To derive the system discussing here, one will need 
more powerful theory. 



2.3 Propagator 

To calculate vacuum energy and graviton mass correction, let us consider concrete forms 
of propagators for small M. In curved spacetime, vacuum is not static or not empty in 
general. We can define that the vacuum is empty at t = by a { p X ' n) \n(0)} = but m a 
later time the background metric mixes the positive and negative frequency components 
[T6] . However, calculations of quantum effects with empty vacuum writing |0) will be 
sufficient because we only try order estimation for latter sections. 
The propagator of empty vacuum can be calculated by 



(2tt) 4 p 2 + m 2 n - (-l)He 
where P livpa = \{ji^r} up + 77^77 

ua ^VpvVo-p) • From now on, we do not explicitly write e. 
As an example, the Feynman rule of the propagator for M = 2 is 

_ z(m 2_ m 2) 

u ^e° (p2 + m g)(p 2 + m?) ^ 1 ] 

This can make the theory renormalizable and the number of courter terms will be finite. 

Since the purpose of this paper is to avoid fine-tuning, M = 2 is not enough. To make 
a super-renormalizable model, we use N = 4 with assuming m 3 = \J wiq — m\ + mf, then 

-i{m\ - ml) (ml - m 2 2 )(m 2 + m\ + 2p 2 ) 
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In the same way, we can get the super- renormalizable model for any even M with the 

condition m^_\ = ^Yln=o(~ l) nm n- Giving one more condition for H > 6, the theory 
becomes finite. For instance, when M = 6, conditions for the finite-field theory are 
m\ — m\ + m\ — m\ + m\ — m\ = and m\{m^ + m 2 ) + m\m\ — m 2 (m 2 + m|) — m\m\. 
In this case, the propagator is 

Dua = + ^ + ^) p (lg) 

\xvpG ^ 2 _|_ m 2)(p2 _|_ m ^(jP _|_ rn 2 )(p 2 + m\)(p 2 + m|)(p 2 + m|) Miypf7 ' 

where /i, /2, /3 denote functions of masses. In the approximation mi « m 3 and m ~ 0, 
it can be written as 

—i(m 2 , — ml) 2 (p 2 + ml)(3(ml — 2ml)p 2 + 2m2 — Smlm 2 ] ) 
(ml — 2m\) 2 (p 2 + ml)(p 2 + m\)(p 2 + m 2 )(p 2 + m|)(p 2 + m 2 i )(p 2 + m 2 ) 

where we used m 2 > m\. With this propagator, all loop calculations become finite and 
we do not need renormalization. 



3 Fermion interaction 

The extension of fermion field to be valid in curved spacetime is formulated via viel- 
bein and covariant derivative. Vielbein can transform the generic metric tensor to the 
Minkowski metric, i.e. = e^e^rjab. The connection for the orthonormal frame is 
written by . Using this connection, the covariant derivative which acts on a spinor 
is given by D^ip = d^ip + \uj ab a a \^ where a ab = — §[7 a ,7 6 ]- We choose the notation of 
gamma matrices to be {7 a ,7 } = —2r] ab . As a general framework, we assume torsion is 
not vanished. Extending the theory based on the Cartan geometry, torsion and curvature 
can be implemented naturally. 

The Lagrangian of fermion field in curved spacetime becomes 



C D = -(-(D^YW + WD^) - m^. (21) 

The connection can be explicitly written as 

1 „,„ „ , 1 „,„ „ , 1 



M ab = ~e/(<9 M e 6i , - d v e hll ) - -e h v (d p e ay - d v e ap ) + -e a p e b a (d a e cp - d p e C(J )e^. (22) 



Since the energy momentum tensor is obtained by T pv = e av T a p = — -^ 8 (£pv^~g) [22] . 
we have 



T>» = - (D^h a e a ^ - -^VD^. (23) 
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This is asymmetric if there is a torsion. Around the Minkowski space background, vielbein 
can be expanded as e a ^ = 5 a ^ — ^h a ^ + • • • and t atl = r\ atl + f /i a/ u + • • • . Then the first 
order of the Hamiltonian density is 

«tet = % ^{d^)l a K\ - ^frfhaWli + ^(d a h° b - d b h° a )^°a ab + a ab ^. (24) 
4 4 1b 

It can make the Feynman rule of 

P =~{E + E')(Yv u ° + tV°) - ^(^V° + ~ P' 3 )^\ (25) 

o 4z 

where S fe = a k <g) 1 2X 2- 





3.1 Newtonian gravity 

In the non-relativistic limit, (0, 0) component of the vertex of fermion-fermion-graviton is 
dominant then 

iK 2 M 1 M 2 P m0Q u\ P l 2 )u{p 2 )u\p' 1 )u{ Pl ) ^ (-1)" 

4 ^-g 2 + m2' W 

where Mi and M 2 are masses of fermions. The calculation of the gravitational potential 
is straightforward: 

v{r) = _^Mi^ { _ ire -^ (27) 

n=0 

Since m is quite light while others are much heavier, it can be approximated by V(r) ~ 
— K MkMi wri ich is 4/3 times larger than Newtonian gravity when G is the observed grav- 
itational constant. This is the famous problem known as the vDVZ discontinuity. 

To see the problem in our framework, let us consider the Einstein equation (JSJ). It 
is sufficient to deal only with hffi since other gravitons do not influence in macroscopic 
scale. The equation of motion of lightest graviton is 

(□ - ml)h^ = - 1 -^{T, V - ilV), (28) 

where we neglect gauge dependent term d^dy/rri^. Considering a point source with T^ u = 
Mr]^Qrj u0 8 3 (x), general solution is 

,, m WnG f d A p 2nMe ip - x 5(p°) , 1 x 



Then 



(0) _ 16nG 2Me~^ (0) _ 16nG Me~ m ° r 
h °° ~ ~ 12^~' " ~k Yh^~ ir [ } 

For the light bending, the angle of deflection can be estimated as 

AGM , . 

a « — , (31) 

where i? is the closest approach to the source. This is the same result of general relativity, 
then we cannot change G to fit the gravitational potential ( )27l) . 

Commonly, it is considered that nonlinear effect of quantum correction can recover the 
smooth connection of massless and massive gravitational theories. Leading contribution 
( )26|) is reliable only when the distance is longer than the Vainshtein radius which is defined 
by (m^Rs) 1 / 5 where Rs is the Schwarzschild radius [23]. This becomes quite long if the 
mass is extremely small, e.g., considering m is the inverse of the Hubble constant, the 
Vainshtein radius of the Sun becomes about lOOkpc which is longer than the size of the 
Milky way. Then nonlinearity is important for the region of the solar system. Summing 
up all the quantum corrections, the result may be correspondent to the case of massless 
graviton (see e.g. [25]). However, this scenario can be applied when corrections from 
loop diagrams are not small compared to the leading order. If usual propagator is used 
and cutoff scale is the Planck mass, loop corrections are comparable. In our case, ghost 
gravitons suppress the corrections so that nonlinear effect is not important. That means 
the problem is serious and the theory is inconsistent with experiments unless we can find 
other remedy. 

As noted in previous section, there may be more fundamental theory that may lead all 
additional terms of the Lagrangian. If the fundamental theory does not involve the prob- 
lem of the vDVZ discontinuity, the theory will be correspondent to the model of massless 
graviton. For instance, basing higher-derivative gravity without Pauli-Fierz mass, the 
discontinuity does not exist [23]. By modifying kinetic terms or strengths of interaction 
terms, we can avoid the problem. A simple example is to take 



£ vDVZ = - -^={h^Uhf v - 2hf u d»d a h^ + 2h$d*ff'hW - h^Dh^)) 

\ 9 (32) 

They can change the equations of motion so that the angle becomes a ~ 16 ^ M . With this 
modification, the gravitational force and light bending are consistent with experiments 
when G = jGn where Gn is the measured gravitational constant. 
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4 Ghost and unitarity 



Ghost particles generally violate unitarity because of negative norms. There are three 
conditions to preserve unitarity [26]: (i) the S- matrix of whole state space V is unitary, 
(ii) physical space V p h ys which is a subset of V is invariant under the S-matrix, namely 
S'Vphys = S^Vphys = Vphys, (hi) physical space has positive semi-definite metric, for any 
|phys) G Vphys, we should hold (phys|phys) > 0. Even when ghost particles are added, all 
of them can be satisfied by assuming the mixing state of normal and ghost particles as a 
fundamental set. Additionally, specific wave packet is assumed to satisfy the unitarity. 

4.1 With renormalizable condition 

We first consider M = 2 so the metric is h uv = h$ + h$. In the analogy of longitudinal 
polarization of photon, we will try to cancel negative norms by constraining the system. 
Suppose the Gupta-Bleuler quantization which gives the constraint (a^ — dp )|phys) = 
where CLp ^ £1X6 annihilation operators of photon, the same constraint for the creation 
operators of gravitons will be useful (ap '^ + a p 1,A ^)|phys) = 0. In this constraint, the 
ghost graviton is always appeared with the normal graviton, then the negative norm will 
not appear. This is the key to satisfy the third condition of unitarity. For the case of 
gauge group, general proof of unitarity is given by [26J and it will be partly applicable for 
the case here. 

The constraint for the physical state will differ from the one of photon due to different 
masses and wave functions. To be more precise, we define 



where wave packets $p are assumed to be J d 3 p $p $p = 1 so that (h u n J\h ( > Tn ^ u ) = 

(—l) n+1 6 nm . Using them, we get the zero-norm state \h nv ) = ^=(|/i]S} + \h$)) where 
the normalization factor is put to be compared with positive-definite norms. The state 
\fipj) can be a physical state if wave packets $ p ^ and $p have the relation J d 3 p(§p^ — 
= 0. Then the wave packets of normal and ghost particles will be the same. The 
constraint is rather simple and seems reasonable if the cancellation works between them, 
nonetheless, the problem is the dependence on the S'-matirx. As it cannot satisfy the 
conditions of unitarity, the constraint needs to be modified. Alternatively, we take (h^u + ^ + 
^A^ + '')|phys) = where 



The new constraint can realize the independence on the S'-matrix since any interaction 




(33) 




(34) 



appears with h$ + h^}. If \h uu ) 



is a physical state, wave packets are related with 
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/ d 3 p($° p /JE { p ] -Q/yjEP) = 0. When the wave packets are Gaussian $J, n) - e ~ P /2 " n 



the relation becomes 



(v^^)3/2 ' 



.1 1 mg, r l 1 



^4'l'4 1 "" ,P| f"fi| 1 = 9 ' (35) 

where U[a, b, z] is the confluent hypergeometric function of the second kind. Assuming 
m /oJo <C 1 and m\/u\ <C 1, we get 

^ 4y^r 2 [|]a; 2 

(Ji « . (36) 

nrrii 

In this way, we can get S\h^ u ) = by the constraint of wave packets. 

Usually, time-independence is supposed as a necessary condition for the conservation 
of unitarity. In fact, if there is a time evolution of \h^ u ), there appears new state \h') = 

\h$) — \h$) and the state of gravitational field turns out to be the mixing of \h^ v ) and 
\h! ). Similar to the neutrino oscillation, this phenomenon occurs since normal and ghost 
gravitons have different masses. Although both states have zero norms, their product 
{hnv\h'^ u ) is negative. Then the appearance of whichever \h^ u ) or \h'^ u ) is problematic. 
Once one can exist, the other is derived by the oscillation. However, it does not a matter 
since the probability to produce \h^ u ) is zero while \h') cannot satisfy the constraint of 
physical state. Therefore the problem can be avoided if we add the constraint that {h^) 
does not exist in the initial condition. 

A proof of unitarity conservation will follow by subtracting zero-norm states from the 
entire physical space if they do not affect for physics. After this, only positive definite 
metrics remain. Nonetheless, a confusion may occur since the energy of the state {h^) is 
not zero. The energy is calculated by 

(KAH\ h n = \ + | , 

where K is the modified Bessel function of the second kind. This energy will be related 
to when the constraint of physical state begins. Below the energy scale, there is no 
problem since \h^}) has positive norm and the ghost \h$) does not appear. In addition, 
\h^u) indicates the ordinary gravitational wave which is indirectly observed. Therefore, 
it is necessary to assume that the physical state of the graviton is first and later 

replaced by \h^ u ) when the energy is beyond the scale that \h^ u ) can appear in the real 
world. The exact value of the lowest energy should be calculated by varying ujq and u\ 
with the constraint of Eq. (135 j) . A subtlety remains since the energy changes with the 
normalization factor which cannot be determined due to zero norm. In other words, the 
zero-norm state affects for the observation indirectly that is completely different from the 
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case of longitudinal photon. For now, as the determination of the factor seems not a 
crucial issue, we just postulate that normalizations of zero norms are the same to those 
of positive definite metrics and then there is no ambiguity. 

The phenomenon of the change of the state from the one of gravitational wave to 
the zero-norm state makes important prediction. Since S-matrix with zero-norm states 
become trivial, all the probabilities including the states are zero. Therefore, the process 
that increases the energy of gravitational wave beyond the energy about mi/2 when the 
state is changed to \h av ) will be canceled out. That implies the energy of gravitational 
wave will have the upper bound. Concretely, there is no way to increase the energy beyond 
the bound. This prediction is not so significant for gravitons, but if we apply the same 
scenario for other elementary particles like electron and neutrino, it is noteworthy. As 
discussed in section 6, particles of the standard model may have ghost partners so that 
the prediction to have upper energy bound will be interesting. 



4.2 With super-renomalizable condition 

We use the same techniques for M > 2 with conditions which can lead higher power of 
propagator. For instance, when M = 4, the super-renormalization condition is m 3 = 

yj TOq — m\ + mf. In this case, the mass m 2 should be heavier than m 3 as m < m\. 

(2) 

Then, we cannot use the cancellation of the negative norm of by constraining to 

the zero norm j5(\h$) + \h<$)). Instead, we require two steps for the physical state: the 

first is the same to M = 2, i.e. (hfy + ^ + /i/^ )|phys) = at the scale larger than the 
minimum of Eq. fl37|) with the state ^(\h$) + \ h^})), the next is to replace this state by 

= 71 En=o \hiMv) from its lowest energy about (mi+m 2 + m 3 )/4 with the constraint 

Y n ^™' + ' ) |phy s ) = 0- To suppress the ghost state a condition for a mass appears, 

precisely, m 3 should be larger than the lowest energy of \h uu ). It can be calculated by 
searching the minimum of 

e m3/2w2 m 2^m| 



n=0 



2 

with varying the parameters u n and taking the condition Yn(~ ^-) nuJ nU[\, —\, ^-] = 0. 
When m 2 and m 3 are similar values, U2 and cj 3 are approximately equal. Then the 
constraints for the masses can be approximated as < m\ < m 2 / y/2 and m 2 < m 3 . 

Let us discuss in more general context. The super-renormalization condition for gen- 
eral M is vrLjsi^x = a/XX~ l)" m n- Ghost masses are assumed to be larger as the number 
is increasing, i.e. m 2n -i < m 2n +i for any n. Defining \h [ ™ ] ) = ^= Yn=o 1 we 

requre that states are changing in a stepwise fasion from \h$) to \hfH), • • • , = \h ul ) 

depending on the energy scale where each state can appear. A constraint for the phys- 
ical state is given by Y^=o ^-M"' + ' > |phy s ) = 0- All ghosts masses should have 
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m 2N+i > 2N+2 En^o" 1 m « to cancel each negative norm by involved in zero-norm states. 
Except norms are zero and unaffected by the S"-matrix so that nothing appears in 

the physical world. 

For individual ghost state, there is at least one condition to be involved in zero-norm 
states, on the other hand, normal gravitons do not need to be a part of such states. For 
instance, the mass of h { $ in M = 6 can be lighter than the lowest energy of \h^ v ). There 
is no other constraint for h$ then \h$) can appear in the final state. The behavior 
of this particle is uncommon because the available momentum is limited, namely, the 
range is from m 4 to | m «- Depending on a parameter region, the motion of h{tv can 
be restricted only in non-relativistic regime. Further, as it is neutral from any gauge 
interaction, life-time will be long enough. Then it is a new candidate of a dark matter. 

When J\[ > 6, all quantum effects become finite by the two conditions: m^-i = 

\/Zn=o(-l) n < and E!:;(-l) m (-El»n + O = 0. Even in this case, the 
above scheme can be applied to preserve unitarity with the constraints hfy +S> |phys) = 
and m 2 N+i > 2 n+2 Sn=o^ m «- The auowe d region of mass parameters will be more 
constrained, but it is possible to obtain the physical state. This will be useful to make 
quantum effects significantly suppressed. 



5 Vacuum energy 

The vacuum energy diverges in the standard calculation of quantum field theory and 
it induces the cosmological constant problem. Using the ghost particle, the propagator 
gains higher power so that the divergence becomes weaker. Using the cutoff scale A, the 
vacuum energy without ghost is of order A 4 and it can be improved by ghost particles 
with the super-renormalization condition into the value of order the fourth power of 
the ghost masses. Experimentally, the cosmological constant is observed about 29meV 4 . 
As we can insert any value for the ghost mass, it is not difficult to lead the observed 
value. However, additional gravitons enforce to change the gravitational potential to 
GivMiMa ^ + v^(_i)n e -nm»r) i The strongest bound of the Newtonian gravity for the 
type QmMiMi^i _|_ e -m y r^ j g roughly m ^ > l/20/im [13]. Then, we cannot use the mass 
scale derived from the cosmological constant A 1 / 4 ~ 1/85/im despite the natural choice 
to be consistent with the observation. Actually, it is not an easy task to find allowed 
parameter region satisfying all the constraints of experiments and unitarity. We abandon 
the complete analysis but show an example of consistent solution. 
The vacuum energy can be calculated by 

^=^<0|GH0>. (39) 

07TG 
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To the leading order, it becomes 



T%= lim V(0|T(^)(x)^^(y))|0>. (40) 

n 

This can be rewritten as 

= hm V I ^ e ^-y) (41) 

When M = 2 and the momentum cutoff is taken, it is integrated into 
T oo = 8K ~ O^ 2 T u = 8(mj - m 2 ) A 2 

128vr 3 ' vac 3 g 47r 3 • ^ ) 

Considering A is the Planck scale, m\ — m 2 , should be O(10 _38 )GeV. Apparently, m\ 
cannot be so small to realize large-scale gravity. On the other hand, if we consider M = 4 
with the condition of super-renormalization, the energy-momentum tensor becomes 

rpfiu = f d*p -i{m\ - ml) (ml - ml) (ml + m 2 2 + 2p 2 )p»p v 
vac J (2tt) 4 (p 2 + m 2 )(p 2 + m 2 )(p 2 + m 2 )(p 2 + m 2 ) ' 1 } 

The leading contribution is proportional to the logarithm of A: 

l07T d 777-2 

It has a wrong sign and the vacuum energy becomes negative. This is because the contri- 
bution from ghost particles is larger than the one of normal gravitons. Although M = 4 
is the least number to be super-renormalizable, it cannot yield a consistent result. 

The simplest predictive model is M = 6 with the super-renormalization condition. To 
get a rough estimation, we take m ~ and m 2 ~ m 4 , then the energy- momentum tensor 
in vacuum becomes 

rpiiv „ _■ f d4 P ~2((m| - m 2 ) 2 - m 2 m 2 )p 4 + 3m 2 777 2 777 2 (p 2 + m|) ^ 
vac ~ 7 (2tt) 4 (p 2 + m 2 )(p 2 + m 2 )(p 2 + m 2 0(p 2 + m 2 ,) ■ °J 

The leading is 

-4-7T 2 A 

T Z « (2^)i(K + ™1 - - ™X) M— ]»T + " " " • (46) 

This time, it is possible to keep a positive energy within the constraint of unitarity. In 
order to get T 00 « 29meV 4 , masses are estimated around 3meV if they have the same 
order without parameter tuning. On the other hand, if we tune the parameters e.g. 
m! ~ 777 3 « 777 2 /a/3, masses can be in any order depending on the level of tuning. More 
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Figure 1: Left figure shows the parameters sets which satisfy the cosmo logical constant 
and the constraints of unitarity. The solid line indicates the mass of m\. Large masses can 
appear when the conditions of finite-field theory are approximately satisfied since they 
suppress the leading order proportional to the logarithm of the cutoff scale which is taken 
as the Planck mass. In right figure, we describe the modified gravity at short distance 
using one parameter set with relatively large masses obtained from the left figure. It is 
compared with the Yukawa-type potential. 



concretely, giving masses around the condition of finite-field theory, the magnitude of 
the leading contribution can be suppressed. For the general case ( l4"5j) . we performed 
numerical analysis to find that all masses can exceed 20meV and Fig. 1 (left) exhibits 
some parameter sets. Taking a set of masses, we predict the modified Newtonian gravity 
by 

V{r) PS GnM\M.2 ^mpr _ e m ir _|_ ^m 2 r _ e m z r _|_ ^n±r _ e m 5 r^ 

r 

One interesting example is (mi, m 2 , m 3 , m 4 , m 5 ) m (21.3, 37.1, 21.4, 30.5, 37.3)meV. Com- 
pared with the Yukawa-type modification, the effective mass is approximately 1/14- 
1/10 jum" 1 which is shown in Fig. 1 (right). Increasing the tuning level, larger mass 
sets will appear whereas they are less interesting since the problem of lightest graviton 
mass arises as discussed in the next section. 

The cosmological constant problem is much improved compared to the fine-tuning 
solution without ghost particles. However, if future experiments verify the inverse square 
law in shorter range, we need to increase the tuning level depending on the range. Strictly 
speaking, as the scale 85/im of natural prediction is already rejected, it does not solve 
the cosmological constant problem at all. Using larger Af, we can take the finite-field 
conditions exactly, but similar tuning will be necessary. It may be important to consider 
additional mechanism or symmetry to reduce the amount of vacuum energy. 
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6 Graviton mass 



There are many models that predict modification of gravity in short range by using 
extra dimension. In our model with ghost particles, it can be distinguished from them 
by predicting modified Newtonian gravity at the cosmological scale via lightest graviton 
mass. From the loop calculations, parameters of gravitons masses are related to the 
lightest mass. The estimation of lightest mass is also important as it is near to the 
experimental bound. Conversely, if the mass of lightest graviton can be measured by 
observation, the parameter set of other gravitons masses can be restricted. After finding 
modified gravity in either short or long range, the model would make stronger prediction. 

The most important effect of loop calculation of graviton mass correction comes from 
the four-gravitons vertex. The Feynman rule of the vertex is given by 



?Ct = —^(h^Kph^doKo + ■■■. (48) 
The Feynman rule of the first term is estimated as 




= -^{--K 2 {p 2 + PsT 1 (p 4 ) V lM V 2At V 3At V 4 ° + 1 2 3 <r> 4). (49) 
P4 4! 2 

Using this rule, the graviton mass correction becomes 

Am 2 „ _ l _ K 2 f dtp 2 y (50) 

The estimation of whole diagrams is difficult and we do not try here. In quantum gravity 
without ghosts, it becomes Am 2 ~ A 2 which needs strong fine-tuning. When M = 2, it is 
Am 2 ~ K 2 A 2 m 2 . This case also needs fine-tuning to make the graviton mass light to be 
consistent with the solar system. To maintain the gravity for the scale larger than lpc, 
Am should be less than 10 _20 meV. 

If a model is super-renormalizable, we can get a realistic graviton mass. The correction 
for M = 4 is 

a 2 _1 2 f -Hm\{m\ - m 2 2 )p 2 pl 

m 2 h J (2tt) 4 {p 2 + m 2 ) (p 2 + m 2 ) {p 2 + m 2 ) (p 2 + m 2 ) ' 1 ' 

Then the propagator becomes 

D -iP^ pu m 2 (m 2 - m 2 )(m 2 + p 2 ) 2 

[jjispcr p2^2 _|_ m 2^(j)2 _|_ m%){p 2 + mjj) + m 2 (mf — m^)^m\ + p 2 )Am 2 

Since Am 2 ~ n 2 m\ ln[A 2 /m 2 ], it is easy to get Am < lpc -1 when ghost masses are 
chosen to fit the cosmological constant. Inserting m\ = A 1/4 and A the Planck mass, it 
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is of order 4 x 10~ 31 meV or 6 x 10 _11 pc _1 . The calculation for the case M = 6 is also 
straightforward. In previous section, one possible solution is given as example. Using 
the parameter set, we found Am ~ 1.4 x 10 _28 meV or 2.2 x 10 _8 pc _1 . Since the case of 
higher-derivative theory is independent from the DGP model (Dvali-Gabadadze-Porrati), 
the graviton mass limit is X g > lOOMpc [15]. Then the solution of the model is marginal 
to the experiments. Taking the parameter set which has larger graviton masses, it will 
be inconsistent with the observation of the cosmological structure. Hence, ghost masses 
are favored to be in the range of 10-100meV to suppress the lightest graviton mass. 
Above calculations reflect just one contribution of the quantum correction and actually 
there are thousands of contributions from similar terms and diagrams. It is inferred that 
by summing up other contributions lightest mass exceeds the experimental bound, then 
we will need some tuning with the bare mass rriQ. The model does not indicate strong 
predictions for modified scales of gravity but implies favored ranges of the modification. 

Lastly, let us comment about divergences induced by interactions with the matter 
sector. They are usually quartic divergences and non-renormalizable. There need infinite 
counter terms to cancel all the divergences even when higher curvature terms are added. 
One possible solution is to consider supergravity [27J and it is known that supergravity can 
reproduce finite results for up to three-loop diagrams. But in general loop, no one knows 
whether the result can be finite. A different approach is to consider higher-derivative 
propagator for all the matter fields then it can make the theory renormalizable with 
power counting arguments. Taking the same procedure of section 4, unitarity can be 
preserved when ghost particles are introduced for any other particle. A critical difference 
is that particles in the standard model are investigated accurately for high energy so that 
ghost masses must be very large. For example, cosmic ray experiments of neutrino find 
that neutrino can have the energy 10 12 GeV, then the mass of ghost partner of neutrino 
must be larger than this value. This mass scale is too high to render the loop correction 
of graviton mass small enough. Although the method is more powerful than supergravity, 
it is still impossible to derive the consistent mass correction of lightest graviton without 
the fine-tuning. There is also a problem when the vacuum energy is calculated for the 
particles of the standard model. If ghost masses substantially exceed the energy scale of 
milli-electron volt, the cosmological constant cannot be explained. Thus, there needs more 
powerful method that can suppress the vacuum energy and the graviton mass correction 
when matter sector is added. For instance, by combining higher-derivative theory and 
supergravity, some strong prediction may be found. In this vein, it there is a strong tool, 
one can predict modified gravity via the quantum effects. 

7 Conclusion 

We have invented new techniques that can preserve unitarity for the system including 
ghost partners. They are so powerful that quantum gravity with matter fields can be 
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renormalizable if ghost partners are introduced for every quantum field. Prediction of 
modified gravity in both short and long ranges is present for the model M = 6 with the 
super-renormalization condition, assuming the contribution of pure gravity is dominant 
and tuning level is weak. It can avoid current experiments covering the inverse square 
law of Newtonian gravity. If future experiments can find modified gravity in short range, 
it can be a clue of ghost partner as well as extra dimension. A difference from the models 
with extra dimension is that it can predict a modification for the gravity on the scale of 
galaxy group and cluster. 

The scheme with ghost partner is indeed powerful, but so far, it is not used widely 
as the unitarity violation is crucial for physics. We have resolved a main difficulty of the 
problem and it can be a good approach to the hierarchy problem. By choosing proper 
physical state with setting the relation with wave packets, positive-semi definite norm 
and independence on the S-matrix are realized. As we did not provide a strict proof 
of the unitarity, there might be some faults. In particular, it is curious that zero-norm 
states indirectly affect for physics by giving the upper bound of energy. This is not seen 
in other quantum theories and possibly meets a problem. All the things we have done 
in this paper are weak statement and prediction of possible new physics toward quantum 
gravity in low energy. 
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